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Sciences  Division,  and  Mr.  Lawrence  A.  Gambino,  Director,  Computer  Sciences 
Laboratory. 

COL  Edward  K.  Wintz,  CE,  was  Commander  and  Director  and  Mr.  Walter  E.  Boge 
was  Technical  Director  of  the  Engineer  Topographic  Laboratories  during  the 
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ANALYSIS  OF  EDGE  DETECTION  ALGORITHMS  ON  DIAL 


INTRODUCTION 


As  part  of  an  effort  to  deliver  to  the  Defense  Mapping  Agency  (DMA)  in 
FY85  a  seraiautomated  feature  extraction  system,  the  Engineer  Topographic 
Laboratories  (ETL)  is  conducting  a  feature  extraction  study.  One  of  the 
objectives  of  this  study  is  to  determine  whether  selected  edge-finding 
algorithms  can  be  used  to  find  such  lineal  features  as  roads  and  boundaries. 


PURPOSE 


The  purpose  of  this  research-  riutH~  is  to  report  on  progress  to  date  on  an 
analysis  of  several  edge-detection  algorithms  that  exist  on  ETL's  Digital 
Image  Analysis  Laboratory ^TdIAL)v image-processing  facility.  In  addition  to 
evaluating  the  algorithms,  this  study  investigated  the  possibilities  of 
integrating  the  algorithms'  resultant  edges  with  other  feature  classification 
algorithms.  Three  fundamental  edge-detection  methods  were  compared: 


1.  A  nonlinear  edge-enhancement  operator. 

2.  A  directional  edge  operator. 

3.  An  extraction  of  the  zero-crossings  from  filtering  an  image  with 

the  Laplacian  of  a  symmetric  Gaussian  function.  _  ■ 


DIAL  FACILITY 


DIAL  is  an  interactive  image-processing  facility  consisting  of  a  host 
computer  (CYBER  170-730),  an  associative  array  processor  (STARAN),  and  a  real¬ 
time  processing  and  display  subsystem.  Several  of  the  edge-detection  program 
modules  (PM's)  require  computationally  bound  mathematical  operations  that 
would  be  too  laborious  and  time-consuming  on  a  sequential  machine  such  as  the 
CYBER  170-730.  In  order  that  the  computations  be  done  as  efficiently  and  as 
expediently  as  possible,  the  services  of  the  array  processor,  STARAN,  which 
processes  data  in  parallel  rather  than  sequentially,  were  engaged.  The  STARAN 
in  conjunction  with  the  CYBER  and  the  DIAL  system  enables  an  operator  to 
produce  interactive  tests  in  a  timely  fashion.  The  DIAL  facility  is  also 
interfaced  with  a  Dunn  camera  system  that  photographs  a  digital  image  on  the 
DIAL  display  screen.  This  camera  was  used  to  take  the  pictures  included  in 
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this  report,  and  although  this  camera  does  not  produce  the  sharpest 
representation  of  the  display  image,  its  photographs  are  an  improvement  over 
any  other  process  we  have  used  thus  far.  A  more  detailed  overview  of  the  DIAL 
facility  is  presented  by  Norvelle. 

TEST  IMAGES 


In  order  to  compare  the  results  of  the  various  algorithms,  all  tests  were 
performed  on  the  same  scenes,  CBP54A  (figure  1)  and  its  stereomate  CBP56A 
(figure  2),  both  of  which  are  subsets  of  a  panchromatic  image  taken  near  the 
U.S. -Canadian  border.  Each  scene  is  a  512  by  512,  8-bit  image. 


Figure  1.  Scene  CBP54A. 


V.  R.  Norvelle,  Interactive  Digital  Correlation  Techniques  for  Automatic 
Compilation  of  Elevation  Data,  U.S.  Army  Engineer  Topographic  Laboratories, 
Fort  Belvoir ,  Va.,  ETL-0272,  October  1981,  AD-A109  145. 
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Figure  2.  Scene  CBP56A  (Stereomate  of  CBP54A) , 


THE  ALGORITHMS 


There  were  three  basic  types  of  edge  operators  tested  in  this  study.  The 
first  was  a  nonlinear  edge-enhancement  operator;  the  second  was  a  directional- 
edge  operator;  and  the  third  employed  a  second-difference  operator  and  then 
determined  the  zero-crossings  of  the  response. 

Nonlinear  -  The  Sobel  operator  was  selected  from  several  nonlinear  edge- 
enhancement  operators  for  comparison  with  other  types  of  operators;  It 
utilizes  a  3  by  3  window  of  information  around  the  pixel  rather  than  only  the 
neighbors  to  the  right  and  below  as  is  the  case  in  a  2  by  2  operator.  The 
Sobel  operator^  computes  a  new  pixel  gray  value  (0)  in  the  following  manner: 


where 


2W.  K.  Pratt,  Digital  Image  Processing,  a  Wi 1 ev-Intersc ience  Publication,  .John 
Wiley  &  Sons,  New  York,  1978. 
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X  =  (A2  +  2A3  +  A4)  -  (Aq  +  2A7  +  Ag) 
and  Y  =  (Ag  +  2A ^  +  A2)  —  (Ag  +  2Ag  +  A4 ) 
and  A^,  i  =  0,  7  are  located  around  the  pixel  P: 


4) 

A1 

a2 

A7 

P 

a3 

*6 

a5 

a4 

Computations  for  the  Sobel  PM  are  performed  by  the  CYBF.R  170-730  computer. 

Directional  Edge  -  The  directional-edge  operator  tested  under  this  study 
is  a  relaxation  technique5  whereby  a  set  of  directional-edge  masks  are 
utiLized  to  generate  a  set  of  likelihood  values  that  an  edge  exists  in  one  of 
eight  directions.  These  Likelihood  values  are  then  enhanced  by  a  relaxation 
process  that  iteratively  updates  the  likelihood  values  in  an  attempt  to 
reliably  define  edge  definitions. 

Because  ot  computational  limitations,  only  edge  masks  of  size  3  by  3  and 
')  by  b  were  used  in  this  study.  The  3  by  3  set  of  masks  in  figure  3  is 
typical  of  directional-edge  masks. 
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Figure  3.  Compass  gradient  Prewitt  edge  masks. 


*k.  a.  Lot l,  R.L.  Faiss,  S.C.  Stadelman,  ST ARAN  Image  Process i ng ,  Final 
’technical  Report,  ■  .S.  Army  Kngineer  Topographic  Laboratories,  Fort  Bel  voir , 
r.i'L  o.l  .  J,  ivtoher  1 9 8 ‘ » ,  AD- At >9  1  928. 
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20.  Zero-crossings  of  7  by  7  Marr  window  with  a 
convolved  with  CBP54A. 


re  21.  Zero-crossings  of  15  by  15  Marr  window  wi 
o=2  convolved  with  CBP54A. 
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Figure  18.  Relaxation  on  CBP54A  using  Kirsch  directions 
masks,  1st  iteration. 


Figure  Id.  Relaxation  on  CBP54A  using  Kirsch  directional 
masks,  3rd  iteration. 


Figure  16.  Relaxation  on  CBPS^A  using  Prewitt  directional 
masks,  1st  iteration. 


Figure  17.  Relaxation  on  CBP51JA  using  Prewitt  directional 
masks,  3rd  iteration. 


When  used  in  the  relaxation  algorithm,  results  of  four  convolutions  (the 
N,  NW,  W,  and  NE  masks)  are  merged  to  form  a  set  of  probabilities  that  are 
"relaxed"  and  iterated  while  using  the  directional  probabilities  of  the 
pixel's  neighbors.  Figure  16  is  the  first  iteration  result  of  the  convolution 
and  relaxation  of  Prewitt  edge  masks  with  CBP5HA,  while  figure  17  shows  the 
third  iteration.  Figures  18  and  19  are  the  first  and  third  iterations  of  the 
same  technique  used  on  the  Kirsch  set  of  directional  masks. 

A  5  by  5  set  of  directional  masks  by  Barnard  was  also  tried  with  the 
relaxation  technique  but  produced  results  no  better  than  the  Kirsch  masks  and 
involved  almost  three  times  the  number  of  computations.  In  order  to  use  the 
results  of  the  relaxation  algorithm  in  a  feature  extraction  scheme,  the 
results  would  again  have  to  be  segmented  into  a  binary  representation  to 
define  an  "edge/no  edge"  existence. 

Figure  21  is  the  result  of  the  Marr  Window  in  figure  7  convolved  with  the  same 
image . 

Figure  22  depicts  the  results  of  a  Marr  Window  with  0=3  convolved  with 
CBP54A. 

In  addition  to  road  and  field  edges,  the  Marr  zero-crossings  algorithm  was 
tested  on  a  drainage  pattern  that  is  a  subscene  of  a  Phoenix  test  model. 
Results  of  this  series  of  tests  are  shown  in  appendix  A. 


DISCUSSION 


This  discussion  will  be  limited  to  the  quality  of  the  results,  not  the 
amount  of  time  elapsed  in  the  computations.  Time  discussions  are  neither 
appropriate  nor  fair  because  some  processes  utilize  a  parallel  processor  for 
computations  while  others  use  sequential  processors.  Then,  too,  the  DIAL 
configuration  utilizes  a  conglomeration  of  computer  hardware,  and  each  edge- 
detection  algorithm  uses  its  own  combination  to  achieve  the  final  results. 


Sobel  -  The  Sobel  operator  performed  admirably  in  delineating  edges 
and  missed  very  few  edges.  It  was  difficult,  however,  to  extract  a  binary 
image  by  thresholding  and  still  keep  the  major  edges  intact.  Clutter  removal 
tended  to  break  up  several  of  these  prominent  edges. 

Directional  Edge  -  Although  the  relaxation  technique  produced  poor 
results  with  the  first  iteration,  the  third  iteration  did  only  slightly  poorer 
than  the  Sobel  operator.  The  Prewitt  set  of  masks  did  not  perform  as  well  as 
the  Kirsch  set  of  masks,  the  latter  giving  finer  detail  along  edges. 


Figure  14.  CBP54A  convolved  with 
Prewitt's  north  gradient  mask. 


Figure  15.  CBP54A  convolved  with  Prewitt’s  west  gradient  mask 
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Figure  13.  Sobel  image  CBSOB  thresholded  at  gray  shade  235. 


Directional  Relaxation  -  As  previously  stated ,  the  relaxation  process 
involves  a  set  of  four  directional  masks.  The  set  included  above  is  „ 

attributed  to  Prewitt."^  Two  other  sets  have  been  studied  here:  the  Kirschr 
directional  masks  and  a  set  of  5  by  5  masks  by  Barnard.  Figure  14  shows  a 
convolution  of  the  Prewitt  north  gradient  mask  with  CBP54A  and  figure  15  is 
the  convolution  of  the  west  gradient  mask  with  the  same  image. 


^H.C.  Andrews,  Interactive  Real-Time  Adaptive  Nonlinear  Convolution,  Image 
Processing  Application  Note,  University  of  Southern  California,  Los  Angeles, 
CA  90007,  no  date. 

°Ibid. 

9S.  Barnard,  Application  of  Image  Understanding  to  Cartography,  SRI 
International  report  8718,  prepared  for  DMAAC/LOP,  March  1981. 
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Figure  12.  Sobel  image  CBSOB  thresholded  at  gray  shade  135. 
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Marr  results  depicted  in  this  report  were  generated  by  a  convolution  of 
the  Marr  windows  with  the  image  in  the  frequency  domain  using  the  FOURIER  PM 
written  by  Mr.  Robert  Pazak  of  ETL.  This  PM  was  modified  to  handle  filters  up 
to  size  49  by  49  on  images  of  size  512  by  512  and  utilizes  the  CYBER. 


RESULTS 


It  is  rather  difficult  to  measure  how  well  an  edge-detection  algorithm 
performs  since  the  results  are  all  visual  and  there  are  no  statistics 
generated  nor  any  standards  with  which  to  compare  the  results.  Too,  the 
laborious  task  of  rectification  and  registration  of  the  images  was  not 
undertaken,  so  it  was  not  possible  to  compare  extracted  edges  to  a  map  sheet. 

Because  we  do  not  intend  to  use  the  results  of  edge-detection  algorithm 
extractions  as  a  final  product,  we  are  not  concerned  by  the  generation  of 
noise  and  clutter.  Those  can  be  removed  through  a  clean-up  process  on  DIAL 
called  RASTER.  In  RASTER,  edges  can  be  thickened  to  join  broken  line 
segments  or  thinned  to  the  center  line  of  line  widths  several  pixels  thick, 
and  clutter  or  noise  can  be  eliminated.  Future  additions  to  RASTER  will 
include  a  directional  line-growing  technique  and  area-fill  for  use  with  the 
Marr  algorithm.  We  are  concerned,  however,  with  how  well  the  algorithm 
depicted  the  presence  of  edges,  if  it  generated  non-existing  edges  or  ignored 
edges  that  do  exist.  Algorithm  results  also  have  to  be  put  into  some  useable 
form  (binary  preferred)  in  order  to  use  them  in  combination  with  other 
feature-extraction  results.  Results  of  the  edge-extraction  algorithms  are 
herewith  presented. 

Sobel.  The  Sobel  operator  yields  an  8-bit  image  representation  upon 
completion.  The  results  are  shown  in  figure  9,  CBS0B1 .  Figure  10,  CBS0B1S, 
shows  the  results  of  the  Sobel  operation  on  the  stereomate  CBP56A.  Note  that 
some  edges,  e.g.,  on  field  in  center  of  photo,  are  not  delineated  as  well  as 
in  figure  9.  Several  lines  of  data  in  the  upper  portion  of  CBS0B1S  are  the 
results  of  scratches  in  the  original  digitized  transparency  and  should  be 
ignored.  The  Sobel  results  were  then  thresholded  in  order  to  utilize  them  in 
conjunction  with  other  feature  extraction  results.  Figures  11,  12,  and  13 
show  figure  9,  a  256  gray  shade  representation,  thresholded  at  gray  shades 
200,  135,  and  235  respectively. 


.  Friend,  Analysis  of  Interactive  Image  Cleansing  via  RASTER-Processing 
Techniques,  U.S.  Army  Engineer  Topographic  Laboratories,  Fort  Belvoir,  Va., 
ETL-0347,  November  1983,  AD-A141  772. 
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Figure  7.  Marr  window  for 


derivative  ol  the  distribution  peaks,  i.e.,  wherever  zero-crossings  in  the 
second  derivative  occur.  The  second  derivative  used  in  our  studv  was  derived 
from  th<_  Gaussian  distribution  in  two  dimensions: 

_r_2_ 

2o2 


">  ■>  2 
where  r~  =  x“~  +  y 


C(r)  =  - 

2  it  a  2 


The  second-order  differential  operator  that  Narr  (and  we)  used  was  the 
Laplacian  because  it  is  orientation-independent,  being  able  to  detect  uniform- 
intensity  changes  running  down  the  v-axis  as  well  as  along  the  x-axis.  "he 
Laplacian  is  defined  as  the  sum  of  the  second  derivatives  of  a  function  with 
respect  to  its  directions: 


3  2g  3  2(; 

LP  (x,y)  =  -  —  + 

3x  2  3v2 


ind  after  the  appropriate  substitutions  are  made: 

G(r) 


LP(x ,y)  = 


X2  +  y  2 


-9 


This  LapLacian  is  used  to  generate  a  filter  window  to  be  convolved  with  the 
image,  using  values  of  x  and  y  to  indicate  the  distance  from  the  center  of  the 
filter.  For  example,  for  a  window  of  size  7,  the  center  filter  element  would 
be  calculated  using  x=y=0,  whereas  filter  element  1,1  would  be  calculated 

using  x=y=-3.  Note  that  since  all  x's  and  y's  in  LP  (x,v)  and  G(r)  are 

squared,  all  filter  windows  generated  are  symmmetric  with  respect  to  both 
axes.  Another  variable  in  the  Laplacian  formula  is  the  spread  parameter, 
which  Marr  utilized  to  deal  with  changes  that  occur  at  different  scales. 

Marr's  theory  in  dealing  with  more  than  one  spread  parameter  was  that  if  two 
or  more  filters  with  varying  spread  parameters  produced  results  that  agreed, 
then  most  likely  the  edges  in  the  image  were  real  (as  owing  to  a  change  in 

reflectance,  illumination,  depth,  or  surface  orientation).  If  the  "edges"  did 

not  occur  in  more  than  one  convolution,  then  most  likely  that  "edge"  was  the 
result  of  noise.  For  a  further  discussion  of  Marr  windowing,  techniques,  see 
appendix  A. 


Figures  f> ,  7,  and  8  are  examples  of  Marr  windows  created  using  several 
spread  parameters. 


In  the  second  step  of  the  relaxation  algorithm,  the  gray  shade  value  of 
each  pixel  is  replaced  by  the  probability  vector  for  that  pixel,  which  is  then 
enhanced  by  a  relaxation  process  using  the  *■  -formation  of  the  probability 
vectors  of  four  neighboring  pixels.  See  figure  5. 


Figure  5.  Neighbors  of  the  reference  pixel. 

A  "best  link"  value  is  computed  between  each  of  the  directional 
probabilities  of  the  reference  pixel  and  all  of  the  directional  probabilities 
of  the  four  neighbors.  The  pair  giving  the  largest  "link"  value  is  used  to 
update  the  directional  probability  of  the  reference  pixel.  According  to  the 
linking  formulation  and  update,  successive  iterations  should  cause  the  true 
edges  to  grow  as  the  number  of  iterations  increases.  Several  methods  of 
display  were  attempted  but  had  drawbacks,  since  it  was  difficult  to  delineate 
direction  on  the  display  screen.  Lott^  settled  on  a  display  of  the  magnitudes 
of  the  largest  probability  measure  of  the  four-component  vector. 

The  relaxation  algorithm  utilizes  the  services  of  the  CYBER  and  the 
STARAN,  the  latter  doing  the  bulk  of  the  computations. 

Zero-Crossings  -  This  algorithm,  developed  by  Marr  and  Hildreth'’,  looks 
for  intensity  changes  in  an  image  by  examining  the  second  derivative  of  the 
normal  distribution.  Intensity  changes  will  occur  wherever  the  first 


A 

R.  W.  Lott,  R.O.  Faiss,  S.G.  Stadelman,  STARAN  Image  Processing,  Final 
Technical  Report,  U.S.  Army  Engineer  Topographic  Laboratories,  Fort  Belvoir, 
Va.,  ETL  0243,  October  1980,  AD-A091  928. 

JD.  Marr  and  F,.  Hildreth,  "Theory  of  Edge  Detection,"  AI  Memo  No.  S18, 
Massachusetts  Institute  of  Technology,  Cambridge,  April  1979. 
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Compass  names  indicate  the  gradient  direction  of  maximum  response.  For 
example,  the  west  mask  produces  a  maximum  output  for  horizontal  reflectance 
changes  from  right  to  left.  The  actual  direction  of  the  line  is  90°  clockwise 
from  the  gradient  direction.  For  a  west  gradient  direction,  the  line 
direction  will  be  to  the  north.  See  figure  4. 


Each  edge  mask  is  convolved  with  the  same  neighborhood  (3  by  3  in  the 
Prewitt  masks)  and  the  response  is  divided  by  the  estimates  of  noise  of  that 
3  by  3  neighborhood  for  normalization  purposes.  The  resultant  value  (one  for 
each  mask)  is  an  estimate  of  the  likelihood  that  an  edge  exists  in  that 
direction.  In  practice,  only  four  gradient  directional  masks  were  used  (NE, 

N,  NW,  and  W) ,  because  a  negative  reflectance  response  indicates  that  the 
gradient  direction  is  actually  in  the  diametric  direction;  e.g.,  a  negative 
response  to  the  west  mask  indicates  an  east  gradient  direction.  At  the  end  of 
the  first  stage  of  this  relaxation  algorithm,  therefore,  there  exists  for  each 
pixel  a  vector  of  four  probabilities  that  an  edge  exists  in  the  -45°,  0°,  45°, 
and  90°  line  direction. 


90°line  direction 


Figure  4.  Relationship  of  line  directions  to  gradient  directions. 


5 


Figure  22.  Zero-crossings  of  27  by  27  Marr  window  with  o=3 

convolved  with  CBP54A. 


Marr  Zero-Crossings  -  The  Marr  algorithm  produces  impressive  visual 
results,  but  it  is  difficult  to  produce  an  "edge/no  edge"  map  from  it.  If  one 
desires  a  "quick  and  dirty"  representation  of  an  edge  map,  then  a  single  Marr 
zero-crossings  representation  in  conjunction  with  the  RASTER  PM  will  produce 
acceptable  results.  The  RASTER  PM,  when  finished,  will  be  capable  of 

eliminating  clutter,  filling  closed  areas,  and  thinning  to  center  lines,  thus 

creating  some  semblance  of  an  edge  map. 

The  smaller  windows  that  are  produced  from  the  smaller  sigma  values 

generate  noisier  images  than  do  the  larger  windows,  but  their  results  present 

a  more  definitive  representation  of  the  edges,  particularly  evidenced  in  the 
samples  depicted  in  appendix  A.  When  two  or  more  windows  with  differing 
sigmas  produce  coincident  edges,  those  edges  most  likely  are  due  to  the 
presence  of  intensity  changes  in  the  image.  Thus,  one  can  produce  an  "edge/no 
edge"  map  by  comparing  several  windows,  but  the  results  will  not  be  as 
complete  as  when  deriving  the  edge  map  from  the  convolved  results  of  a  small 
window. 
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CONCLUSIONS 


1.  The  nonlinear  method  of  edge  detection  generates  a  good  visual 
representation.  It  is  difficult  to  use  the  results,  though,  in  combination 
with  other  feature  extractors,  as  pertinent  data  is  lost  converting  the 
results  to  a  useable  format. 

2.  The  images  produced  from  the  directional  edge-relaxation  method  of 
detection  were  inferior  to  those  produced  by  the  Sobel  method.  Edges  were  not 
as  well  defined  and  the  images  were  of  poorer  quality. 

3.  The  zero-crossings  operator  produces  a  unique  representation  of  edges, 
but  its  utility  will  depend  upon  developing  an  effective  cleanup  PM  like 
RASTER. 

4.  Use  of  a  stereo  image  is  valuable  in  picking  up  additional  information 
because  of  the  different  perspective  presented  by  the  stereomate.  However,  in 
order  to  transfer  newly  identified  data  from  one  image  to  its  stereomate,  it 
is  necessary  to  register  the  two  images  in  order  to  extract  and/or  overlay  the 
additional  information. 


APPENDIX  A 


RESULTS  OF  VARIOUS  MARR  WINDOW  CONVOLUTIONS 


Background  -  The  Marr  algorithm  was  originally  tested  at  ETL  on  a  drainage 
pattern  subset  of  a  Phoenix  model  and  the  results  were  impressive.  However, 
the  Marr  windows  produced  by  the  larger  sigma  values  (o  >_  3)  were  quite  large 
and  required  an  excessive  amount  of  computer  time  to  perform  the  convolution 
of  the  Marr  window  with  the  image.  In  an  effort  to  decrease  the  turnaround 
time,  it  was  decided  to  perform  the  convolutions  in  frequency  space  instead  of 
in  the  spatial  domain  as  was  originally  done.  This  step  successfully  reduced 
computation  time  from  30  minutes  to  70-90  seconds  for  larger  sized  filters. 

It  was  also  decided  to  examine  the  filter  sizes  to  determine  if  the  pertinent 
information  could  be  contained  in  a  reduced  filter  size,  and  what  filter  sizes 
are  appropriate  to  the  various  sigma  values.  This  appendix  presents  the 
results  of  the  convolutions  on  the  Phoenix  model  with  varying  Marr  filters. 

Theoretically,  the  sum  of  the  filter  values  should  be  equal  to  zero.  But 
according  to  the  table  in  figure  A1 ,  for  a  sigma  value  of  3,  any  of  the 
filters  from  size  23  by  23  up  to  size  33  by  33  could  qualify,  depending  on  the 
accuracy  required.  From  the  table  labeled  figure  A2,  it  would  appear  that  the 
outer  five  rows  would  contribute  insignificantly  to  the  results  because  their 
values  are  relatively  small.  The  sum  of  the  filter  "weights"  for  this  filter 
is  equal  to  .0001. 

Results  -  It  can  be  seen  from  the  table  in  figure  A1  that  by  extracting  a 
smaller  window,  e.g.,  a  19  by  19,  from  the  center  of  the  entire  Marr  filter, 
the  convoluted  result  becomes  more  positive,  thus  biasing  the  results  to  some 
extent.  This  bias  is  due  to  the  fact  that  the  eliminated  filter  values  are 
all  negative.  The  effect,  though,  can  be  an  advantage  as  can  be  seen  in 
figures  A4  through  A17,  which  are  the  results  of  convolutions  of  Marr  windows 
with  the  drainage  pattern  subset  of  our  PHOENIX  model,  figure  A3- 
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a 

WINDOW 

FILTER  RUM 

1 

5 

.934 

1 

7 

.051 

1 

9 

.001 

2 

1 1 

9.524 

2 
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Figure  A1 .  Sum  of  filter  values  for  various 
Marr  windows. 
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Figure  All.  Zero-crossing  of  21  by  21  Marr  window, 

o=3. 
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Figure  A17.  Zero-crossing  of  31  by  31  Marr  window, 

a  =  5. 


After  looking  at  the  results  in  figures  A5,  A7 ,  and  A12  through  A15,  one 
sees  that  little  can  be  done  with  extracting  edges  even  by  means  of  clutter 
removal.  Please  note  that  for  these  filters  the  sum  of  the  weights  approaches 
zero.  In  the  series  of  pictures  for  o  =  3,  a  better  choice  for  a 
representative  pattern  of  edges  would  be  figure  A10,  where  a  cleanup  process 
might  be  more  manageable.  Note  that  for  figure  A10,  the  sum  of  the  filter 
values  is  16.513. 

Additional  exercises  were  undertaken  to  force  the  sum  of  the  filter  values 
to  be  zero  when  the  filter  was  large  and  when  a  smaller  filter  with  which  to 
work  was  desired.  Two  options  were  tried:  (1)  by  subtracting  from  each  value 
of  the  filter  an  average  correction  based  on  the  sum  of  the  filter  weights, 
and  (2)  by  decreasing  only  the  positive  values  of  the  original  filter. 

Neither  option  was  acceptable,  as  the  resulting  zero-crossings  images  were 
just  as  "busy"  as  the  original  and  did  not  suit  our  purposes  for  extracting  an 
edge.  Another  exercise  involved  a  thresholding  technique.  The  floating  point 
array  that  is  created  by  the  convolution  was  histogrammed  in  order  to  get  some 
idea  of  the  distribution  of  the  values.  From  this  histogram,  a  threshold 
range  could  be  selected  to  test  the  existence  of  a  zero-crossing  in  order  to 
make  the  procedure  less  sensitive  to  noise  in  the  image.  It  was  found, 
however,  that  thresholding  destroyed  the  well-defined  closed  ellipses, 
rendering  the  image  useless  for  some  rastering  techniques. 
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